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Abstract 

Let u be a function of n independent variables x , ...,x n , and U = (uij) the 
Hessian matrix of u. The symplectic Monge- Ampere equation is defined as a linear 
relation among all possible minors of U. Particular examples include the equation 
det U = 1 governing improper afflne spheres and the so-called heavenly equation, 
^13^24— 1*231*14 = I, describing self-dual Ricci-flat 4-manifolds. In this paper we clas- 
sify integrable symplectic Monge-Ampere equations in four dimensions (for n = 3 
the integrability of such equations is known to be equivalent to their linearisability) . 
This problem can be reformulated geometrically as the classification of 'maximally 
singular' hyperplane sections of the Pliicker embedding of the Lagrangian Grassman- 
nian. We formulate a conjecture that any integrable equation of the form F(v,ij) = 
in more than three dimensions is necessarily of the symplectic Monge-Ampere type. 
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1 Introduction 



Let us consider a function u(x l , ...x n ) of n independent variables and introduce the n x n 
Hessian matrix U = of its second order partial derivatives. The symplectic Monge- 
Ampere equation is a PDE of the form 

M n + M n _ 1 + ... + M 1 + M = (1) 

where Mi is a constant-coefficient linear combination of all I x I minors of the matrix 
U, < I < n. Thus, M n = detU = Hessw, M is a constant, etc. Equivalently, these 
PDEs can be obtained by equating to zero a constant-coefficient n-form in the 2n variables 
x\ Ui. Equations of this type belong to the class of completely exceptional Monge- Ampere 
equations introduced in [6]. Geometric and algebraic aspects of symplectic Monge- Ampere 
equations have been thoroughly investigated in [23, 3]. We point out that the left hand side 
of (1), M(U), constitutes the general form of null Lagrangian densities, that is, functionals 
of the form J M{U)dit which generate trivial Euler-Lagrange equations [2]. The class of 
equations (1) is invariant under the natural contact action of the symplectic group Sp(2n), 
which is thus the equivalence group of our problem. All subsequent classification results 
will be formulated modulo this Sp(2n)-equivalence. 

In the case n = 2 one arrives at the standard Monge- Ampere equations, 

e(u u u 22 - u\ 2 ) + au u + f3u l2 + ^u 22 + 5 = 0; (2) 

these are known to be the only equations of the form F(un, Ui 2 , u 22 ) = which are 
linearisable by a transformation from the equivalence group Sp(4). 

The case n = 3 is also understood completely: one can show that any non-degenerate 
symplectic Monge- Ampere equation is either linearizable, or Sp(6)-equivalent to either of 
the canonical forms, 

Hess u = 1, Hess u = Uu + u 22 + m 33 , Hess u = Un + u 22 — m 33 , (3) 

see [23, 3]; we point out that all three canonical forms are Sp(6)-equivalent over C. The 
first equation arises in the theory of improper affine spheres, while the second describes 
special Lagrangian 3-folds in C 3 [8, 19]. Here the non- degeneracy is understood as follows: 
let F(uij) = be a symplectic Monge-Ampere equation. Consider the linearized equation, 
dF/duij Vij = 0, obtained by setting u — > u + ev and keeping terms of the order e. The 
non- degeneracy means that the corresponding symbol, dF/du^ defines an irreducible 
quadratic form. 

The problem of integrability of symplectic Monge-Ampere equations was addressed in 
[15] based on the method of hydrodynamic reductions [12, 13, 14]. Without going into 
technical details of this method, let us formulate the main result needed for our purposes: 

Theorem 1 [15]. A non- degenerate three-dimensional symplectic Monge-Ampere equa- 
tion is integrable by the method of hydrodynamic reductions if and only if it is linearisable. 
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In particular, the PDEs (3) are not integrable. Although this result is essentially negative, 
it will be crucial for the classification of integrable equations in higher dimensions (where 
the situation is far more interesting). Before we proceed to the description of the main 
results, let us clarify the geometry behind symplectic Monge- Ampere equations (1) and 
the linearisability/integr ability conditions. Let us consider the Lagrangian Grassmannian 
A, which can be (locally) identified with the space of n x n symmetric matrices U. The 
minors of U define the Pliicker embedding of A into projective space P N (we identify A 
with the image of its projective embedding). Thus, symplectic Monge- Ampere equations 
correspond to hyperplane sections of A. For n = 3 we have A 6 C P 13 , and linearisable 
equations correspond to hyperp lanes which are tangential to A 6 . Therefore, for n = 3 the 
linearisability condition coincides with the equation of the dual variety of A 6 , which is 
known to be a hypersurface of degree four in (p 13 )* (we refer to [21, 22, 26] for a general 
theory behind this example). In Sect. 2 we provide the following characterisation of 
linearisable equations in any dimension: 

Theorem 2 For a non- degenerate symplectic Monge- Amp ere equation (1) the following 
conditions are equivalent: 

(1) The equation is linearisable by a transformation from Sp(2n) . 

(2) The equation is invariant under an n 2 -dimensional subalgebra of Sp(2n). 

(3) The equation corresponds to a hyperplane which contains an osculating subspace O n _2 
of the Lagrangian Grassmannian of the order n — 2. 

For n = 3 the third condition reduces to the requirement of tangency. 

In Sect. 3 we address the problem of integrability of symplectic Monge- Ampere equa- 
tions in four dimensions, n — 4. Among the best known four-dimensional integrable 
examples one should primarily mention the 'heavenly' equation, 



[29] which is descriptive of Ricci-flat self-dual 4- manifolds [1]. It was demonstrated in [13] 
that this equation is integrable by the method of hydrodynamic reductions. Although, 
in principle, the method of hydrodynamic reductions can be applied in any dimension, it 
leads to a quite complicated analysis. One way to bypass lengthy calculations is based 
on the following simple idea: suppose our aim is the classification of four-dimensional 
integrable equations of the form (1) for a function ^(a: 1 , x 2 , x 3 , x A ). Let us look for traveling 
wave solutions in the form 



where Q is an arbitrary quadratic form in the variables x l , x 2 , x 3 , x A . The substitution 
of this ansatz into (1) leads to a three-dimensional symplectic Monge- Ampere equation 
which must be integrable for any values of constants a,/3,j, and an arbitrary quadratic 



- ^23^14 = 1 



(4) 




or, more generally, 



u = u^x 1 + ax 4 , x 2 + (3x 4 , x 3 + 7a; 4 ) + Q(x, x), 
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form Q. Since, in three dimensions, the integrability conditions are explicitly known (and 
are equivalent to the linearisability) , this provides strong restrictions on the original four- 
dimensional equation which are therefore necessary for the integrability. In fact, in the 
present context they turn out to be sufficient: if all three-dimensional equations obtained 
from a given four- dimensional PDE by traveling wave reductions are linearisable, then 
the PDE is integrable. The philosophy of this approach is well familiar from the soliton 
theory: symmetry reductions of integrable systems must be themselves integrable. Thus, 
for the first heavenly equation, traveling wave solutions are governed by 



which is a three-dimensional symplectic Monge- Ampere equation. One can show that it is 
indeed linearisable for any values of constants. This approach suggests a simple geometric 
characterisation of integrable equations in four dimensions. Let us first point out that for 
n = 4 the Lagrangian Grassmannian A 10 C P 41 is foliated by a 7-parameter family of A 6 
where each A 6 corresponds to a collection of Lagrangian planes in the symplectic space 
V 8 which pass through a fixed vector. 

Theorem 3 A non-degenerate four-dimensional symplectic Monge-Ampere equation is 
integrable if and only if the corresponding hyperplane is tangential to the Lagrangian 
Grassmannian A 10 along a ^-dimensional subvariety which meets all A 6 C A 10 . 

Thus, integrable equations correspond to maximally singular hyperplane sections. The 
classification of such hyperplanes leads to a complete list of integrable equations in four 
dimensions: 

Theorem 4 Over the field of complex numbers any integrable non- degenerate symplectic 
Monge-Ampere equation is Sp(8)- equivalent to one of the following normal forms: 

1. Mn — U22 — W33 — W44 = (linear wave equation); 

2. W13 + W24 + — u±2 — (second heavenly equation); 

3. M13 = ^12^44 — ^14^24 (modified heavenly equation); 
4- u i3 u 24 — U14U23 = 1 (first heavenly equation); 

5. Mn + U22 + M13M24 — Wi 4 m 23 = (Husain equation); 

6. awi 2 M34 + (3ui 3 ii24 + 7^14^23 = (general heavenly equation), a + j3 + 7 = 0. 

Eqs. 2-6 are non- linearisable, and contact non-equivalent. They have appeared in different 
contexts in [29, 18, 23]. To the best of our knowledge, the integrability of Eqs. 3 and 6 was 
not recorded before. We refer to [34, 35, 33, 16, 10, 11, 27, 31, 4, 5] for the Hamiltonian, 
twistorial and symmetry aspects of heavenly-type equations. 

Symplectic Monge-Ampere equations play a distiguished role in the classification of 
all integrable PDEs of the 'dispersionless Hirota' type, 



a(«i2«i3 - U11W23) + /3(wi3«22 - M12M23) = 1, 




(5) 
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which constitute a single relation among second order partial derivatives of a function 
u{x l , x n ). For n = 3 integrable equations of this type were studied in [15, 28, 32], 
revealing a remarkable correspondence with the theory of hypersurfaces of the Lagrangian 
Grassmannian, generalised hypergeometric functions and GL(2, Restructures. Among 
them, only linearisable equations belong to the class (1). On the contrary, for n > 4 all 
known integrable examples belong to the class (1). This suggest the following conjecture: 

Conjecture. For n > 4 any non- degenerate integrable equation of the form (5) is neces- 
sarily of the symplectic Monge-Ampere type. 

Provided this conjecture is true, Theorem 4 gives a complete list of integrable equations 
of the form (5) in four dimensions. 



2 Characterisation of linearisable equations: proof of 
Theorem 2 

The results of this section are valid for an arbitrary number of independent variables 
n. We begin with a few general remarks. The Lagrangian Grassmannian A is a variety 
of dimension n(n + l)/2 which can be locally identified with n x n symmetric matrices 
U = (uij) so that Uij can be viewed as local coordinates on A. The action of Sp(2n) is 
defined by the infinitesimal operators 



Xi 







duij ' 



d d 



OUjc OU; 



_ d \- d 



ou ss du sk 

s s=f=k 

notice that = Xji and Pij = Pji, while ^ Lji. Thus, we have n{n+ l)/2 operators 
Xij, n 2 operators and n(n + l)/2 operators P^. Altogether, they generate a Lie algebra 
Sp(2n) of dimension n(2n + 1). We will work with the affine Pliicker embedding of the 
Lagrangian Grassmannian A in the projective space P N ~ l specified by the position vector 

r = (Mi,M 2 ,..,M n ); 

here Mi denotes a collection of all Z x Z minors of U, and N = C 2 n n - C 2 n + 2 is the total 
number of minors of U (equivalently, N is the dimension of the vector space of effective 
n- forms). Thus, in the case n = 3 one has 

dimMi = 6, dimM 2 = 6, dimM 3 = 1, A 6 C P 13 . 

Similarly, for n = 4 one has 

dimMx = 10, dimM 2 = 20, dimM 3 = 10, dimM 4 = 1, A 10 C P 41 . 
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With each point of the Pliicker embedding of A one can associate a sequence of osculating 
subspaces 

1 c 2 c 3 c ... c O n = P N ~ 1 

where 0\ is the tangent space, O2 is the second osculating subspace, etc. At the point 
corresponding to = the subspace 0\ is spanned by vectors from M 1; the subspace 
O2 is spanned by vectors from Mi U M 2 , etc. A symplectic Monge- Ampere equation is 
naturally identified with a hyperplane section of A. Below we present three equivalent 
characterisations of linearisable equations. 

Theorem 2. For a non- degenerate symplectic Monge-Ampere equation in n dimensions 
the following conditions are equivalent: 

(1) The equation is linearisable by a transformation from Sp(2n). 

(2) The equation is invariant under an n 2 -dimensional subalgebra of Sp(2n). 

(3) The equation corresponds to a hyperplane which contains an osculating space O n -2- 

Proof: 

equivalence (1) (2): since all n- dimensional linear equations are equivalent under 
the (complexified) action of Sp(2n), let us consider the Laplace equation, 



Mil H \-Unn = 0. 



A direct calculation shows that it is invariant under the n 2 -dimensional subalgebra of 
Sp{2n) generated by the operators (i 7^ j), Xu—Xjj, L^—L^ (i 7^ j), and La. Since 
this property is manifestly Sp(2ri)-invariant, any non-degenerate linearisable equation is 
invariant under a subalgebra of Sp(2n) of dimension n 2 . The converse is true for any (not 
necessarily Monge-Ampere) non-degenerate equation of the form F(iiij) = 0. Indeed, 
let G be a symmetry group of such equation. We can always assume that the point 0, 
specified by ity = 0, belongs to the hypersurface in the Lagrangian Grassmannian A 
corresponding to our equation. Let G be the stabilizer of this point in G. Note that 
dimG — dim Go < dim A — 1 = n(n + l)/2 — 1, as G takes the equation to itself. The 
stabilizer of the point G A under the action of Sp(2n) has the form: 



P 



A AB 




A e GL(n),B = l B 



Since the equation is non-degenerate, we can always bring it to the canonical form: 

F(uij) = ttu H \-u nn + o{uij) = 0. (6) 

This form (together with the point 0) is stabilized by the elements of P given by the 
condition A 6 CO(n), B = 0. Hence, dim Go < dim CO(n) so that dimG < dim A — 1 + 
dim CO{n) = n 2 . The equality holds only if Go = CO{n). But the subgroup of scalar 
matrices in CO(n) acts by non-trivial scalings on the terms of the order 2 and higher in 
(6). Hence, if Go = GO(n), all higher order terms vanish identically, and our equation 
coincides with the Laplace equation. 
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equivalence (1) -<=>- (3): Let us fix a point on A corresponding to the parameter values 
Uij = (since Sp(2n) acts transitively on A, all points are equivalent). The osculating 
space O n -2 at this point is specified by the equations M n _i = 0, M n = 0. Any hyperplane 
containing this osculating space corresponds to an equation which is a linear combination 
of minors of U of the orders n — 1 and n. Any such equation linearises under the Legendre 
transformation. This finishes the proof. 

Remark. Linearisable equations of Monge- Ampere type are of interest in their own, for 
instance, the equation 

u tt (l + u xx + u yy ) - u 2 xt - u 2 yt = e 
governs the evolution of 'Kahler potentials' [9]. It corresponds to the Lagrangian 

J [u 2 (l + u xx + Uyy) + 2eu\ dxdydt, 

and linearises under the partial Legendre transformation, x = x, y = y, t = u t , u — 
tu t — u, u x = —u x , Uy = —u y , ui = t, taking the form 

1 "ii ^yy ^^ii' 

Notice that all symplectic Monge-Ampere equations are Lagrangian. This is a general 
fact which does not require the linearisability/integr ability. 

3 Classification of integrable equations 

In this section we classify integrable four-dimensional symplectic Monge-Ampere equa- 
tions up to the action of the equivalence group Sp(8). In Sect. 3.1 we begin with a ge- 
ometric discussion which establishes a correspondence between integrable equations and 
'maximally singular' hyperplane sections of the Lagrangian Grassmannian A 10 C P . 
This leads to the main classification result proved in Sect. 3.2. Lax pairs, symmetry alge- 
bras and singular varieties associated with integrable equations arising in the classification 
are discussed in Sect. 3.3 - 3.5. 

3.1 Geometric preliminaries: proof of Theorem 3 

Let v G V be a vector in the 8-dimensional linear symplectic space V. Lagrangian 
subspaces passing through v form a sub-Grassmannian A 6 C A 10 C P 41 . Thus, the 
Lagrangian Grassmannian A 10 is foliated by a 7-parameter family of A 6 . Each A 6 lies in 
some P 13 C P 41 . Given a four- dimensional symplectic Monge-Ampere equation (1), let 7r 
be the corresponding hyperplane in P . The condition that all traveling wave reductions 
of the four- dimensional equation are linearisable is equivalent to the requirement that the 
corresponding hyperplane it is tangential to all A 6 C A 10 . This can only happen if either 
of the following conditions is met: 

(1) The 7-parameter family of sub-Grassmannians A 6 possesses an envelope E C A 10 
which is contained in 7r; 
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(2) The hyperplane n is tangential to A 10 along a 4-dimensional subvariety X 4 C A 10 
which meets each A 6 (we point out that no hyperplane can be tangential to A 10 along a 
subvariety of more than four dimensions). 

The first case can be ruled out straightaway: the 7-parameter family of sub-Grassmannians 
A 6 is 'too big' to possess a non-trivial envelope (e.g., a one-parameter family of lines in 
the plane always possesses an envelope, however, the full two-parameter family of lines 
doesn't). Thus, we are left with the second possibility, which establishes the necessity part 
of Theorem 3. The sufficiency will follow from the classification of all non- degenerate equa- 
tions with a four-dimensional variety of tangency provided in Sect. 3.2: the requirement 
of the existence of such X A proves to be very strong indeed, and leads to a finite list of 
examples all of which turn out to be integrable. 

In the following we will call the variety of tangency X 4 the singular variety of the 
symplectic Monge-Ampere equation defined by the hyperplane section tt. Indeed, X A 
coincides with the set of singular points of the intersection of ir with A 10 . 

The condition that the singular variety X A meets each A 6 is essential for the integra- 
bility. In more detail, consider the incidence variety Y C A 10 x PV, 

Y = {(l,(v))\vel}, 

along with the two projections TTi : Y — > A 10 and 7r 2 : Y — > PV. Then ir^ (X A ) is a 7- 
dimensional closed algebraic variety, and the projection 112(^1 1 (X 4 )) is a closed algebraic 
subvariety in PV. The condition that X A meets all sub-Grassmannians A 6 is equivalent 
to 7T2(7t 1 ' 1 (A a4 )) = PV. This allows us to check this condition using affine charts on A 10 
which have non-empty intersections with X A . 

Example 1. Let us describe the singular variety x A for the non-degenerate linear equation 
U12 — W34. As shown above, the corresponding hyperplane has a contact of order 2 with 
the Lagrangian Grassmannian A 10 at the Lagrangian subspace given by x l = 0. There 
is also a natural conformal structure on 1^ associated with this equation. Explicitly, it is 
given by du\du 2 — du^du^ = 0. A direct computation in various affine charts of A 10 shows 
that the set of points where the corresponding hyperplane is tangential to A 10 consists of 
all Lagrangian subspaces / such that one of the two following cases holds: 

• I has a 3-dimensional intersection with l^; 

• I has a 2-dimensional intersection with which is isotropic with respect to the 
above conformal structure. 

Altogether, these points form a 4-dimensional closed algebraic variety with three irre- 
ducible components: one for subspaces with a 3-dimensional intersection, and two for 
subspaces with 2-dimensional isotropic intersections. 

It is easy to see that the first component meets all A 6 , that is for any v G V there 
exists a Lagrangian subspace I C V such that v G I and dim / R Zoo = 3. Moreover, if v 
does not lie in l^, such Lagrangian subspace I is unique. 

Example 2. Consider the equation Hessu = 1, which is known to be non-integrable. Let 
us show that its singular variery X A misses most of the A 6 . Partial Legendre transform 
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brings it to 14111*22 — uf 2 = U33U44 — u\ A . The set of points where the corresponding 
hyperplane is tangential to A 10 is given by the equations un = w 22 = W12 = and 
U33 = M44 = M34 = 0. In other words, the intersection Xq of X 4 with an affine chart of A 10 
defined by subspaces x l = and Uj = is represented as the following set of symmetric 
matrices, 

/ !i 13 u 14 \ 
W 23 W 24 
U13 U 23 
\Ul4 U 24 / 

It is easy to check that ^(^^(Xq)) lies in a quadric in PV and, hence, ^(^^(Xq)) 7^ 
PV. This confirms that the equation Hessw = 1 is not integrable, even though the 
corresponding hyperplane is tangential to A 10 along a 4-dimensional variety. 

As we have seen in Example 1, the variety of tangency X 4 may consist of several 
irreducible components with different geometric properties. However, if the equation 
under study is non-degenerate and non-linearisable, any two generic Lagrangian subspaces 
from X 4 must have trivial intersections. This follows from the two Lemmas which analyse 
the cases when generic pairs of Lagrangian subspaces of X 4 have two- or one-dimensional 
intersections. 

Lemma 1. Let X 4 be a four-dimensional subvariety of A 10 such that any two generic 
Lagrangian subspaces from X 4 have two-dimensional intersections. Then there exists a 
fixed Lagrangian subspace L which has three-dimensional intersections with all subspaces 
from X 4 . This situation corresponds to linearisable equations. 

Proof: 

Let li, I2, 13 be any three generic subspaces from X 4 such that dim(Zj D lj) = 2. Then 
dim(/i fl I2 fl h) = 1, and we can introduce the Lagrangian subspace / = span{Zi PI h, h H 
h,h H I3}. Any other generic Lagrangian subspace from X 4 will have one- dimensional 
intersections with Zi Pi Z 2 , I3 and I2 H £3. Thus, it will have a three-dimensional inter- 
section with I. Any hyperplane 7r which is tangential to A 10 along X 4 will automatically 
contain the osculating subspace to A 10 at I, so that the corresponding equation will be 
linearisable. 

Lemma 2. There exists no non-degenerate symplectic Monge- Ampere equation for which 
any two generic Lagrangian subspaces in the singular variety X 4 have one- dimensional 
intersections. 

Proof: 

Let l\ and Z 2 be two Lagrangian subspaces from X 4 such that dim(7i fl I2) = 1. Let 
v be a vector in the intersection. Consider the Lagrangian Grassmannian A 6 C P 13 
generated by v. Then the travelling wave reduction in the direction of v produces a three- 
dimensional linearisable Monge- Ampere equation such that the corresponding hyperplane 
in P 13 is tangential to A 6 at two distinct points, l\ = h/(v) and Z 2 = h/(v)- The 
points l\ and Z 2 correspond to three-dimensional Lagrangian subspaces in V 6 which do 
not intersect by construction. This situation is, however, not possible for linearisable 
equations in three dimensions (even degenerate). Thus, the travelling wave reduction in 
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the direction of v must produce the 'zero' equation, that is, A 6 C 7r fl A 10 where it is the 
hyperplane corresponding to our Monge-Ampere equation. When l\ and I2 are allowed 
to vary within X 4 , one obtains at least a four-parameter variety of pairwise intersections 
v (as explained below). Thus, we obtain a four-parameter family of A 6 which belong to 
the nine-dimensional hyperplane section n fl A 10 . This, however, is not possible since the 
union of any four-parameter family of A 6 must be essentially ten-dimensional. 

It remains to prove that we have at least a four- dimensional variety of pairwise inter- 
sections v . Let us fix l\ and consider the set of all intersections Yil\) = {h n/2, h £ AT 4 }. 
Let s be the dimension of Y(li). Fixing v in Y(li) we obtain a (4 — s)-parametric family 
of Lagrangian subspaces I in X A which contain v. Each of them contains an s-parameter 
subvariety Y(l). Taking a union of Y(l) over all I passing through v, one obtains a 
four-parameter family of pairwise intersections v. This finishes the proof of Lemma 2. 

Let us now return to the general case when any two generic Lagrangian subspaces 
of X 4 have trivial intersections. Take any two such planes, say l\ and I2, and introduce 
canonical coordinates x\ Ui in V such that I1J2 become the coordinate planes x l = 
and Ui = 0, respectively. Suppose now that a symplectic Monge-Ampere equation 
corresponds to a hyperplane it which is tangential to the Lagrangian Grassmannian A 10 
at the points corresponding to l\ and l^- In coordinates x l ,Ui, any such equation takes a 
'purely quadratic' form generated by a linear combination of 2 x 2 minors of the matrix 
U only. Notice that such form is invariant under GL(4) C Sp(8). Let us choose a third 
Lagrangian plane £3 which does not intersect l\ and 1%, and bring it to a normal form using 
the remaining GL(4)-freedom. The requirement that it is tangential to A 10 at the point 
corresponding to I3 imposes further constraints on the quadratic part of the equation. 
Further details of the classification are provided in the next Section. 

3.2 Classification results: proof of Theorem 4 

Theorem 4. Over the field of complex numbers any integrable non- degenerate symplectic 
Monge-Ampere equation is Sp(8)- equivalent to one of the following normal forms: 

1. Mn — U22 — W33 — W44 = (linear wave equation); 

2. Ui 3 + U2A + U11U22 — U12 — (second heavenly equation); 

3. u\s = M12M44 — M14M24 (modified heavenly equation); 
4- U13U24 — U14U23 = 1 (first heavenly equation). 

5. Uu + U22 + ^13^24 — ^14^23 = (Husain equation); 

6. aui2U 34: + f3ui 3 U24 + 7M14U23 = (general heavenly equation), a + (3 + 7 = 0. 

Proof: 

As shown above, we can always assume that the Monge-Ampere equation contains only 
quadratic terms, or, in other words, the corresponding hyperplane n C P 41 is tangential 
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to the Lagrangian Grassmannian A 10 at the points l\ — oo (specified by x l = 0) and l 2 = 
(specified by Ui = 0). Any other Lagrangian plane Z 3 which has trivial intersections with 

11 and Z2 is defined by a non-degenerate symmetric 4x4 matrix U\j. It can be viewed up 
to the natural action of the subgroup GX(4) C Sp(8) which stabilizes l\ and l 2 . Over the 
field of complex numbers all such symmetric matrices are equivalent to each other. Fix l 3 
such that M14 = 1, u 23 = —1 and all other w^- equal 0. Assume that, in addition to l\ and 
l 2 , our hyperplane 7r is also tangential to A 10 at the point Z 3 . 

The set of all quadratic Monge- Ampere equations such that the corresponding hyper- 
planes are tangential to the Lagrangian Grassmannian A 10 at the points h,l 2 ,l 3 forms a 
10-dimensional vector space spanned by the following expressions: 

E = u n u 2 2 - uf 2 , 

Ei = l/2(%% - u 12 u 14 + u 22 u 13 - u 12 u 23 ), 

E 2 = l/6(u n w 44 - u\ A + u 22 u 33 - u 2 23 ) + 1/3(2mi 3 m 2 4 - ^14^23 - ^12^34), 
E 3 = l/2(u 33 u 24 - u 23 u 3A + u 44 u 13 - u u uu), 
E 4 = u 33 u M - U3 4 ; 

and 

F = u u u 33 - u 2 l3 , 

F 1 = l/2(u n u 3 4 - u r3 u u + U 33 Ui 2 - ni 3 n 2 3), 

F 2 = l/6(w n «4 4 - «?4 + u 22 u 33 - ul 3 ) + 1/3(2mi 2 m 3 4 - M14M23 - M13M24), 
F 3 = l/2(u 22 u 34 - u 23 u 2A + u A4 ui 2 - U14U24), 
F 4 = u 22 u 44 - ul 4 . 

Denote by E and F the subspaces spanned by Ei and Fj, respectively. Note that per- 
mutations of indices (1,2) and (3,4) leave both subspaces invariant, while permutations 
(1,4) and (2,3) interchange them. Both subspaces are also invariant under the natural 
action of SO(4,C) given as the set of all linear symplectic transformations stabilizing li, 

1 2 and / 3 . Moreover, it is possible to show that under the identification of 5*0(4, C) with 
(SX(2,C) x SL(2,C))/(±E 2 ,±E 2 ), each of the copies of SX(2,C) preserves the decom- 
position E © F, acting irreducibly on one summand and trivially on the other. 

The standard model for an irreducible 5-dimensional SL(2, C)-action is given as the 
following action on degree 4 polynomials: 

(: s) @3) ■ 

The vectors Ei and Fj are chosen in such a way that the correspondence Ei -H- t l and 
Fi •<->■ t\ i — 0, . . . , 4, establishes the isomorphisms of SX(2, (C)-actions on subspaces E 
and F with the standard model. Polynomials of degree 4 are easily classified with respect 
to the SL(2, C) action and the multiplication by a non-zero scalar. We have the following 
set of canonical representatives according to the roots of the polynomial: 

(0) 0; 
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(1) 1 (all four roots of p(t) coincide and are equal to oo); 

(2) t (p(t) has one triple root oo and one simple root 0); 

(3) t 2 (p(t) has two double roots and oo); 

(4) t 2 — 1 (p(t) has one double root oo and two simple roots ±1); 

(5) (t 2 — l)(t — a), a ±1 (p(t) has four simple roots oo, ±1, a). 

This classification gives us an easy way to classify vectors in E@F viewed up to the action 
of SO (A, C), permutation of E and F and the multiplication by a non-zero scalar. Namely, 
decompose any non-zero vector v in E © F into the difference v e — Vf,v e £E,Vf£F and 
match each of these two vectors with the corresponding polynomial in the standard model. 
This gives a pair of polynomials (p(t),q(t)), which fully describe the vector v. Next, we 
classify this pair up to independent actions of SL(2, C) on each of them, permutations 
p -H- q and a multiplication by a non-zero constant. We easily get the set of normal forms 
v i,j — {Piif)i c 1jif)), i,j = 0,...,5,i>j,c^0, where Pi and qj are polynomials listed 
above under the i-th and j-th items. We can assume since the case c = is included 
in the subcase q(t) = qo(t) = 0. Next, we consider each of these normal forms case-by-case 
and check for which of the corresponding Monge-Ampere equations the variety of points 
where the corresponding hyperplane tt is tangential to the Lagrangian Grassmannian A 10 
has dimension 4. The complete list of all such cases is given below: 



1. 


pit) 


= Q(t) = it 2 - 


l)(f 


2. 


pit) 


= q(t) =t 2 - 


i; 


3. 


p{t) 


= t 2 -l, q(t) 


= t 2 


1. 


pit) 


= q(t) = t 2 - 




5. 


pit) 


= q(t) = t; 




6. 


pit) 


= t, q{t) = l; 




7. 


pit) 


= q(t) = l; 




8. 


pit) 


= t 3 - t, q{t) 


= 0; 


9. 


pit) 


= t, q(t) = 0; 




10. 


pit) 


= 1, q(t) = 0. 





The cases 4, 7 and 10 correspond to the degenerate equations M13M24 — ^12^34 = 0, U\\U2% — 
u i2 = u u u 33 ~ u i3 an d ^11^22 — u i2 = 0) respectively. In fact, the last two equations are 
equivalent to the degenerate linear equations U22 = ^33 and U22 = 0, respectively, via the 
partial Legendre transform of (x 1 ,-^): (x 1 ,-^) H- (— UijX 1 ). 
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The case 9 is given by 



Wll«24 - Ui 2 U U + n 22 '«13 - "12^23 = 0, 

and is equivalent to a non-degenerate linear equation (and, hence, to the linear wave 
equation) via the partial Legendre transforms of {x x ,ui) and (x 2 ,w 2 ). 

Similarly, the case 8 is equivalent to Hessw = 1. Indeed, the polynomial t 3 — t is 
equivalent to t 4 — 1 under the action of SL(2, C) since in both cases all 4 roots are 
distinct and their cross-ratio is equal to —1. The pair of polynomials p(t) = t A — 1, 
q{t) = corresponds to the equation Unu 22 — u\ 2 = ^33^44 — u\ 4 which can be brought to 
Hess u = 1 via the partial Legendre transform of (x 1 , u\) and (x 2 , w 2 ). As we have seen in 
Example 2, this equation is not integrable. 

Other cases can be considered in a similar manner. It appears that they correspond 
to the remaining non-linear integrable Monge-Ampere equations listed in Theorem 4. 
Namely, equation 5 can be brought to the modified heavenly equation via the change of 
variables (x 1 , x 2 , x 3 , x 4 ) 1— > (x 1 , x 2 + x 3 , x 2 — x 3 , x 4 ), partial Legendre transform of (x 1 , U\) 
and a permutation of indices. 

Equations from cases 2 and 3 are brought to the Husain and the first heavenly equa- 
tions, respectively, via the partial Legendre transform of (x l ,Ui). 

The case 6 is brought to the second heavenly equation via the partial Legendre trans- 
forms of (x^iii) and (x 2 ,u 2 ), and a permutation of indices. 

Finally, to establish the correspondence between equation 1 and the general heavenly 
equation one needs to transform the canonical form (t 2 — l)(t — a) to t A + 2bt 2 + 1 under 
the action of SL(2, C) and then use the following change of variables: ( 
(x 1 + x 4 , x 2 + x 3 , x 2 — x 3 , x 1 — x 4 ). 

The obtained correspondence is summarized in the following table: 



Case 


p(t) 




Equation 


1 


(t 2 -l)(t-a) 


(t 2 -l)(t-a) 


general heavenly 


2 


t 2 -l 


t 2 -l 


Husain equation 


3 


t 2 -l 


t 2 


first heavenly 


4 


t 2 


t 2 


degenerate equation 


5 


t 2 


t 


modified heavenly 


6 


t 


1 


second heavenly 


7 


1 


1 


degenerate equation 


8 


t(t 2 - 1) 





Hess u = 1 (non-integrable) 


9 


t 





linear wave equation 


10 


1 





degenerate equation 



It remains to point out that all equations appearing in Theorem 4 are indeed integrable: 
they possess the required number of hydrodynamic reductions (for the first and second 
heavenly equations this has been demonstrated in [13, 14]), along with the Lax pairs 
depending on an auxiliary spectral parameter. This completes the proof of Theorem 4. 

Remark. Although all equations listed in Theorem 4 are not 5p(8)-equivalent (in fact, 
not even contact equivalent), there may exist more complicated connections among them. 
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Below we describe the construction of [29] which links the first and second heavenly 
equations. Let us take the second heavenly equation in the form 9 x i + 6* 22 +^11^22 — = 0> 
here 9 is a function of the four independent variables X . X . X . X . Introducing the two- 
form 

Q = (dx 1 — 9^dx 1 + 6^ 2 <ix 2 ) A (dx 2 + 9^dx 1 — 9ndx 2 ), 

one can verify the relations dfl = 0, Q A Q = 0. Thus, by Darboux's theorem, there exist 
coordinates x 3 ,x 4 such that Q = dx 3 A dx 4 '. This implies the expansions 

dx 1 — 9^dx 1 + 9i2-dx 2 = u^dx 3 + u^dx 4 , dx 2 + ^cfe 1 — 9ndx 2 = U2 3 dx 3 + U2^dx 4 

where Uij are so far arbitrary coefficients satisfying the relation M13M24 — M14M23 = 1- 
Rewriting these relations in the form 

dx 1 = 922,dx 1 — 9i^dx 2 + u Vi dx 3 + u^dx 4 , dx 2 = —9i^dx 1 + 9i\dx 2 + U2 3 dx 3 + U2idx 4 , 

one can see that, by virtue of the relation dx 1 /dx 2 = dx 2 /dx 1 , there exists a function 
^(x 1 , x 2 , x 3 , x 4 ) such that x 1 = du/dx 1 , x 2 = du/dx 2 . In this case coefficients w^- become 
second order partial derivatives of u, and the relation Ui 3 u 24i — M14M23 = 1 becomes the sec- 
ond heavenly equation. The passage from X ■ X - X - X to the new independent variables 
X/ - X/ - oc • rX C9jii be viewed as a mult i- dimensional analogue of reciprocal transformations 
familiar from the (1 + l)-dimensional theory. This link is highly nonlocal, and in any case 
not a contact equivalence. The only 'local' relations one can write here are the following: 

x 1 = ui, x 2 = u 2 , 9-22 = un, 9 U = -u 12 , 9 Vl = u 22 . 

One can expect similar non-local links between other equations listed in Theorem 4. 

3.3 Lax pairs 

In this section we present Lax pairs for integrable equations from Theorem 4. Some 
of them readily follow from the Lax pair for the six-dimensional version of the second 
heavenly equation [34, 30], 

U15 + U 2 6 + U 13 U 2 4 ~ U U U 23 = 0, 

which is given by a pair of vectors fields 

X x = <9 6 + Mi 3 <9 4 - u u d 3 + Xdx, X 2 = d 5 - w 23 <9 4 + u 2i d 3 - A<9 2 , 

A =const, such that [X 1; X 2 ] = modulo the equation. Straightforward dimensional 
reductions provide Lax pairs for the first heavenly, second heavenly and Husain equations 
(we point out that the general heavenly equation does not travelling; wave 

reduction of this six-dimensional equation). 
Second heavenly equation Ux3 + w 2 4 + UxxU22 — u\ 2 = ; 

Xx = d 4 + u u d 2 - Ux2dx + A<9i, X 2 = d 3 - u 12 <9 2 + u 22 <9i - A<9 2 . 
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Modified heavenly equation u 13 = u^u^ — U\ 4 u 24 . 

X 1 = u u d 2 - u 12 d 4 + Xdx, X 2 = -d 3 + u 44 d 2 - u 24 d 4 + A<9 4 . 

First heavenly equation u 13 u 24 — U\ 4 u 23 = 1: 

Xi = u 13 d 4 - u u d 3 + Xdx, X 2 = -u 23 d A + u 2i d 3 - \d 2 . 

Husain equation Un + u 22 + u l3 u 24 — W14W23 = 0: 

= <9 2 + ^13^4 - ^14^3 + A<9i, X 2 = di - u 23 <9 4 + u 2A d 3 - A<9 2 . 

General heavenly equation aui 2 u 3 4 + fiu\ 3 u 2A + 7M14M23 = 0: 

Xi = w 34 <9i - Mi 3 <9 4 + 7A(m 34 <9i - Mi4<9 3 ), X 2 = u 23 d 4 - u 34 d 2 + (3\(u 34 d 2 - u 24 d 3 ). 

Notice that in the last case the commutator of X\,X 2 is not identically zero on the solu- 
tions, this condition is substituted by [Xi,X 2 ] = mod X±,X 2 - Thus, the 2-dimensional 
distribution spanned by X±, X 2 is integrable. Modifications of the inverse scattering trans- 
form and the <9-dressing method for Lax pairs of this type were developed in [24, 25, 5]. 

3.4 Symmetry algebras 

In this section we describe symmetry algebras for equations from Theorem 4. To be 
precise, we consider only those symmetries which belong to the equivalence group Sp(8) 
or, equivalently, stabilisers of the corresponding effective forms. It should be emphasized 
that the full contact symmetry algebras of these equations are infinite dimensional. We 
adopt the notation of Sect. 2. 



Equation 


Generators of the symmetry algebra 


Dimension 


linear 

wave 


X11+X22, X11+X33, X11+X44, X12, X13, X\4, X23, X24, X34, 
L11+L22+L33+L44, L12+L21, L13+.L31, L14+L41, L23 — L32, L24 — L42, L34 — L43 


16 


second 
heavenly 


Xn, X12, X13— X24, X14, X22, X23, X33 — L24, ^34+2I/23i -^44 — L13, 
L12—L43, L21—L34, L14—L23, 2L11+L22+L44, L11+2L22+L33 


14 


modified 
heavenly 


X\\, X22, X23, X24 — L?,4, X33, X34, X44+L32, 
ill, ^22+^33, 2L33+L44, L24, Pll, P44— 21/23 


13 


first 

heavenly 


Xll, X12, X22, X33, X34, X44, L12, 
L21, L34, L43, L11— L22, L33—L44, L11+L22—L33—L44 


13 


Husain 


X\\—X22i X12, X33, X34, X44, Lll+I/22, 
^12—^21, 1-33—^44, £34, ^43, P1I—P22, P\2 


12 


general 
heavenly 


Xii, Pa, i — 1,...,4. 


12 



15 



The table shows that equations from Theorem 4 are not equivalent under the action 
of the symplectic group. Indeed, the second heavenly equation is the only one with a 
14- dimensional symmetry algebra. Both the modified heavenly and the first heavenly 
equations have 13-dimensional symmetry algebras, but only the symmetry algebra of the 
modified heavenly equation has an invariant 2-dimensional non-isotropic subspace in C 8 . 
Similarly, both the general heavenly and the Husain equations have 12-dimensional sym- 
metry algebras, but the symmetry algebra of the general heavenly equation is reductive, 
while that of the Husain equation is not. 

Remark. For any effective n-form u on the 2n-dimensional symplectic space V with the 
symplectic 2-form Q we can define a bilinear form B^ [23], 



B w : (X,Y) ^ 



jjcw A iyUJ A Q 



which will be symmetric for odd, and skew-symmetric for even values of n. It is known 
that is proportional to Q for all even values of n (we thank Bertrand Banos for pointing 
this out). As oj is only defined up to a scalar multiple, the form B^ also makes sense only 
up to multiplication by a non-zero constant. Hence, if B^ = XQ, we can always assume 
that either A = (i.e., the form B^ vanishes identically), or A = 1 and B^ = Q. It is clear 
that the vanishing of A is an invariant condition for the form u. Below we list the values 
of A for integrable Monge- Ampere equations from Theorem 4: 



Equation 


Coefficient A 


linear wave 





second heavenly 





modified heavenly 





first heavenly 


1 


Husain 


1 


general heavenly 


1 



For instance, this confirms that the modified heavenly equation is not contact equivalent 
to the first heavenly equation, although both have stabilizers of dimension 13. 



3.5 Geometry of singular varieties 

We recall that any integrable equation corresponds to a hyperplane which is tangential 
to the Plucker embedding of the Lagrangian Grassmannian A 10 along a four-dimensional 
subvariety X 4 which meets all A 6 C A 10 . The second condition means that the correspond- 
ing congruence of Lagrangian subspaces fills the symplectic space V^ 8 . In this section we 
provide invariant descriptions of A 4 for all examples arising in the classification. 

Linear equations: Fix a Lagrangian subspace L C V 8 . Then A 4 consists of all La- 
grangian subspaces which have three-dimensional intersections with L. We emphasize 
that this constitutes only one irreducible component of the singular variety. See Exam- 
ple 1 for more details. 
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Second heavenly equation: Fix a Lagrangian subspace L and a two-dimensional 
isotropic plane I such that / C L C V 8 . Let k be the P 1 of lines in /, and K the P l 
of three-dimensional subspaces in L which contain I. Fix a projective isomorphism be- 
tween these two P x 's, k ^ K. Let c be the union, over k, of two-dimensional isotropic 
planes which pass through k and are contained in the corresponding K. Similarly, let C 
be the union, over K, of all Lagrangian subspaces which pass through K. Fix a projective 
isomorphism c -B- C which covers the isomorphism k -H- K. Then X 4 consists of all La- 
grangian subspaces which contain an isotropic plane from c, and have a three-dimensional 
intersection with the corresponding Lagrangian subspace C. 

Modified heavenly equation: Fix a decomposition V 8 = V 2 + V 6 into a direct sum 
of two skew-orthogonal symplectic spaces. Let I be a two-dimensional isotropic plane in 
V 6 , and L the corresponding co-isotropic four-dimensional subspace, I C L C V 6 . Let k 
be the P 1 of lines in I, and K the P l of three-dimensional subspaces in L which contain 
/. Fix a projective isomorphism between these two P lJ s, k <f-» K. Finally, let F be the 
family of Lagrangian subspaces in V 6 which pass through one of the fc's, and have a two- 
dimensional intersection with the corresponding K. Notice that F is three-dimensional. 
Then X 4 consists of all Lagrangian subspaces which have one- dimensional intersections 
with V 2 , and intersect V® along a Lagrangian subspace from F. 

First heavenly equation: Fix a decomposition V 8 = V 4 + V 4 into a direct sum of 
two skew-orthogonal symplectic spaces. Fix a two-dimensional Lagrangian plane in each 
copy of V 4 , denote these planes l\ and respectively. Let Fi be the family of all 
Lagrangian planes in the first copy of V 4 which have one-dimensional intersections with 
l\. Similarly, let F2 be the family of all Lagrangian planes in the second copy of V 4 which 
have one-dimensional intersections with I2. Then X 4 consists of all Lagrangian subspaces 
which intersect each copy of V 4 along two-dimensional Lagrangian planes from F\ and 
F 2 , respectively. 

Husain equation: Fix a decomposition V 8 = V 2 + V 2 + V 4 into a direct sum of three 
skew-orthogonal symplectic spaces. Fix a two-dimensional Lagrangian plane / C V 4 and 
consider a family F of all Lagrangian planes in V 4 which have one-dimensional intersec- 
tions with /. Then X 4 consists of all Lagrangian subspaces which have one-dimensional 
intersections with each copy of V 2 , and intersect V 4 along a Lagrangian plane from F. 
General heavenly equation: Fix a decomposition V 8 = V 2 + V 2 + V 2 + V 2 into a 
direct sum of four skew-orthogonal symplectic spaces. Then X 4 consists of all Lagrangian 
subspaces which have one- dimensional intersections with each copy of V 2 (four-secant 
subspaces) . 

One can show that congruences of Lagrangian subspaces in V 8 corresponding to X 4 
are of class one, that is, a generic vector in V 8 lies in a unique Lagrangian subspace from 
X 4 . 

Acknowledgements 

We have benefited from discussions with Bertrand Banos, Maciej Dunajski, James Grant, 
Valentin Lychagin, Joseph Landsberg, Maxim Pavlov, Volodya Rubtsov and Artur Sergeev. 



17 



We thank the LMS for their financial support of BD to Loughborough making this collab- 
oration possible. We also thank the ICMS for giving us an opportunity to spend two weeks 
in Edinburgh under the 'research in groups' program where this work has been completed 
and the investigation of higher dimensional Monge- Ampere equations was initiated. 



References 

[1] M.F. Atiyah, N.J. Hitchin and I.M. Singer, Self-duality in four-dimensional Rieman- 
nian geometry, Proc. Roy. Soc. London Ser. A 362, no. 1711 (1978) 425-461. 

[2] J.M. Ball, J.C. Currie and P.J. Olver, Null Lagrangians, weak continuity, and varia- 
tional problems of arbitrary order. J. Funct. Anal. 41, no. 2 (1981) 135-174. 

[3] B. Banos, On symplectic classification of effective 3- forms and Monge- Ampere equa- 
tions, Diff. Geom. Appl. 19 (2003) 147-166. 

[4] L.V. Bogdanov and B.G. Konopelchenko, On the heavenly equation hierarchy and 
its reductions, J. Phys. A 39, no. 38 (2006) 11793-11802. 

[5] L.V. Bogdanov and B.G. Konopelchenko, On the <9-dressing method applicable to 
heavenly equation. Phys. Lett. A 345, no. 1-3 (2005) 137-143. 

[6] G. Boillat, Sur l'equation generale de Monge-Ampere a plusieurs variables, [On the 
general Monge-Ampre equation in several dimensions] C. R. Acad. Sci. Paris Ser. I 
Math. 313, no. 11 (1991) 805-808. 

[7] CP. Boyer and J.D. Finley, Killing vectors in self-dual Euclidean Einstein spaces, J. 
Math. Phys. 23 (1982) 1126-1130. 

[8] E. Calabi, Improper affine hyperspheres of convex type and a generalization of a 
theorem by K. Jorgens, Michigan Math. J. 5 (1958) 105-126. 

[9] S.K. Donaldson, Nahm's equations and free-boundary problems, arXiv:0709.0184. 

[10] M. Dunajski and L.J. Mason, Twistor theory of hyper-Kahler metrics with hidden 
symmetries. Integrability, topological solitons and beyond. J. Math. Phys. 44, no. 8 
(2003) 3430-3454. 

[11] M. Dunajski and L.J. Mason, Hyper-Kahler hierarchies and their twistor theory. 
Comm. Math. Phys. 213, no. 3 (2000) 641-672. 

[12] E.V. Ferapontov and K.R. Khusnutdinova, On the integrability of (2+l)-dimensional 
quasilinear systems, Comm. Math. Phys. 248 (2004) 187-206. 

[13] E.V. Ferapontov and M.V. Pavlov, Hydrodynamic reductions of the heavenly equa- 
tion, Class. Quantum Grav. 20 (2003) 2429-2441. 



18 



[14] E.V. Ferapontov and K.R. Khusnutdinova, Hydrodynamic reductions of multi- 
dimensional dispersionless PDEs: the test for integrability, J. Math. Phys. 45 (2004) 
2365-2377. 

[15] E.V. Ferapontov, L. Hadjikos and K.R. Khusnutdinova, Integrable equations of the 
dispersionless Hirota type and hypersurfaces in the Lagrangian Grassmannian, In- 
ternational Mathematics Research Notices 2009; doi: 10.1093/imrn/rnpl34; arXiv: 
0705.1774 (2007). 

[16] J.D.E. Grant, On self-dual gravity, Phys. Rev. D 47 (1993) 2606-2612. 

[17] J. Harris, Algebraic geometry: a first course, Springer- Verlag, New York (1992), 
328pp. 

[18] V. Husain, Self-dual gravity as a two-dimensional theory and conservation laws, Clas- 
sical Quantum Gravity 11, no. 4 (1994) 927-937. 

[19] D. Joyce, Lectures on special Lagrangian geometry, in: Global theory of minimal 
surfaces, Clay Math. Proa, 2, Amer. Math. Soc, Providence, RI (2005) 667-695. 

[20] A. Kushner, V.V. Vychagin and V.N. Rubtsov, Contact geometry and nonlinear 
differential equations, Encyclopedia of Mathematics and its Applications 101, Cam- 
bridge University Press, 2007. 

[21] J.M. Landsberg and L. Manivel, The projective geometry of Freudental's magic 
square, Journal of Algebra 239 (2001) 477-512. 

[22] J.M. Landsberg and L. Manivel, Construction and classification of complex simple 
Lie algebras via projective geometry, Sel. Math., New ser. 8 (2002) 137-159. 

[23] V.V. Lychagin, V.N. Rubtsov and I.V. Chekalov, A classification of Monge-Ampere 
equations, Ann. Sci. Ecole Norm. Sup. 26, no. 4 (1993) 281-308. 

[24] S.V. Manakov and P.M. Santini, Inverse scattering problem for vector fields and the 
Cauchy problem for the heavenly equation, Phys. Lett. A 359, no 6 (2006) 613-619. 

[25] S.V. Manakov and P.M. Santini, On the solutions of the second heavenly and Pavlov 
equations, arXiv:0812.3323. 

[26] S. Mukai, Simple Lie algebras and Legendre variety, Preprint, translation of the 
article in Nagoya Suri Forum, 3 (1996) 1-12. 

[27] F. Neyzi, Y. Nutku and M.B. Sheftel, Multi-Hamiltonian structure of Plebanski's 
second heavenly equation, J. Phys. A 38, no. 39 (2005) 8473-8485. 

[28] A.V. Odesskii and V.V. Sokolov, Integrable pseudopotentials related to generalized 
hypergeometric functions, arXiv:0803.0086. 



19 



[29] J.F. Plebanski, Some solutions of complex Einstein equations, J. Math. Phys. 16 
(1975) 2395-2402. 

[30] J.F. Plebanski and M. Przanowski, The Lagrangian for a self-dual gravitational field 
as a limit of the SDYM Lagrangian, Phys. lett. A 212 (1996) 22-28. 

[31] M.B. Sheftel and A. A. Malykh, On classification of second-order PDEs possessing 
partner symmetries, J. Phys. A 42, no. 39 (2009) 395202, 20 pp. 

[32] A.D. Smith, Integrability of Second-Order Partial Differential Equations and the 
Geometry of GL(2)-Structures, PhD thesis, Department of Mathematics, Duke Uni- 
versity, 2009. 

[33] I.A.B. Strachan, The symmetry structure of the anti-self-dual Einstein hierarchy, J. 
Math. Phys. 36, no. 7 (1995) 3566-3573. 

[34] K. Takasaki, An infinite number of hidden variables in hyper-Kahler metrics, J. Math. 
Phys. 30, no. 7 (1989) 1515-1521. 

[35] K. Takasaki, Symmetries of hyper-Kahler (or Poisson gauge field) hierarchy, J. Math. 
Phys. 31, no. 8 (1990) 1877-1888. 



20 



